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Abstract 

This paper addresses alternative procedures to the ESO supphed pipehne pro- 
cedures for the reduction of UVES spectra of two quasar spectra to determine the 
value of the fundamental constant /i = Mp/M^. at early times in the universe. 
The procedures utilize intermediate product images and spectra produced by the 
pipeline with alternative wavelength calibration and spectrum addition methods. 
Spectroscopic studies that require extreme wavelength precision need customized 
wavelength calibration procedures beyond that usually supplied by the standard 
data reduction pipelines. An example of such studies is the measurement of the 
values of the fundamental constants at early times in the universe. This article de- 
scribes a wavelength calibration procedure for the UV- Visual Echelle Spectrometer 
on the Very Large Telescope, however, it can be extended to other spectrometers as 
well. The procedure described here provides relative wavelength precision of better 



* Based on observations made with ESO Telescopes at the La Silla or Paranal 
Observatories under program IDs 68.A-0106 and 70.A-0017 
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than 3 x 10 for the loiig-sht Thorium- Argon cahbration lamp exposures. The gain 
in precision over the pipeline wavelength calibration is almost entirely due to a more 
exclusive selection of Th/Ar calibration lines. 
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1 Introduction 



This paper serves two main purposes. The first is to establish a detailed record 
of the UVES wavelength calibration and spectrum addition procedures used in 



the m easurement of the fundamental constant fi = Mp/M^ ([Thompson et al. 



20081 ) based on spectra taken with UVES on the VLT. The second purpose is 
a description of an alternative to the pipeline method of establishing a wave- 
length cahbration of Very Large Telescope (VLT) UV-visual Echelle Spec- 
trometer (UVES) spectra. In a belief that a scientific result on the value of a 
fundamental constant in the early universe should be described well enough 
that it can be repeated by other researchers, we have provided more detail 
than is usually present in the description of a data reduction method. 

The standard data reduction pipeline provided by the European Southern Ob- 
servatory (ESO) for UVES does an excellent job of providing well calibrated 
spectra for use in most scientific investigations. Certain investigations, how- 
ever, req uire a more precise w avelength calibration than is provided by the 



pipehne. iMurphy et al.l (120071 ) point out that the UVES pipeline wavelength 
calibration has systematic deviations in accuracy that can adversely affect 
precise measurements. In addition, proper validation of the accuracy of re- 
sults obtained from a pipeline data reduction requires an exact knowledge of 
the pipeline data reduction procedures. This is often not easily available or 
discernible for a standard pipeline. One example of such observations is the 
measurement of the values of fundamental constants in the early universe such 
as a, the fine structure constant and fi, the ratio of the proton to electron mass, 
through precise spectroscopy of distant objects. The faintness of such objects 
often makes the use of techniques such as iodine absorption cells impractical. 
This article describes a set of observing and data reduction procedures that 
can improve upon the wavelength accuracy relative to that available with the 
standard UVES pipeline. The procedures achieve a wavelength accuracy of 
AA/A ~ 2.5 X 10"'' for Thorium-Argon lines in the long-slit calibration spec- 
tra for UVES with a slit width of 0.8 arc seconds. The main component in 
achieving this accuracy is the proper selection of appropriate calibration lines. 



Email addresses: rthoinpson@as.arizona.edu, jbechtold@as.arizona.edu, 
John.Black@chalmers.se, C.J.A.P.Martins@danitp.cam.ac.uk (C.J. A. P. 
Martins). 
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similar to that of iMurphy et al.l (120071 ). 



The examples used in this analysis are VLT UVES spectra of Q0347-383 and 
Q0405-443 taken in January of 2002 and 2003 respectivelj0- These particular 
spectr a have been used to d etermine the value of fj, at early tirn e s in t he uni- 
verse ( Reinhold et al. ( 2006 ) and references therein. iKing et aL ( 20081 )). The 
measurements of n are based on the absorption line wavelengths of molecular 
hydrogen produced in damped Lyman alpha clouds along the line of sight to 
the QSO. The redshifts of t he lines are z=3.024 9007 for Q0347-383 and 



z=2.5947361 for Q0405-443 ([Thompson et all . 120081 ). In this work, however 



we concentrate on the long-slit calibration line spectra that accompany these 
observations and the proper addition of the individual 2x2 binned pixel object 
spectra that are an intermediate product of the UVES pipeline data reduction. 
In all of the following the word pixel refers to the 2x2 binned pixels produced 
during readout and maintained during the pipeline processing. The goal is to 
find, as accurately as possible, the true wavelength solution for each order of 
the spectrum and how that solution varies with time during the course of the 
observations. All solutions are for vacuum wavelengths. The wavelength cal- 
ibration techniques are independent of the objects being observed. The best 
results occur when a rigorous sequence of object and calibration observations 
is performed. This sequence was not strictly followed in this case but the 
sequence is close enough to be useful. 



2 The observations 



The observations of Q0347-383 and Q0405-443 with UVES on VLT occurred 
during the nights of January 7-9 2002 f or Q0347-383 and Janu ary 4-6 2003 for 



Q0405 -443 from a program described by lPetitjean et al.l (l2004f ) and llvanchik et al 



( 120051 ). On each of the nights three separate spectra of the QSO were taken 
with accompanying short and long-slit calibration lamp integrations at the 
same grating setting, all with a slit width of 0.8 arc seconds. We note that 
the long-slit calibration lamp spectra are included in the observation defini- 
tion, presumably to prevent a resetting of the grating position between the 
object and calibration spectra. In this analysis we only use the long slit cal- 
ibration spectra which have a higher signal-to-noise ratio than the short slit 
spectra. For Q0347-383 two long-slit calibration spectra were taken each night 
immediately after the first and second object spectra with the third object 
spectrum occurring immediately after the second long-slit calibration spec- 
trum. For Q0405-443 three long-slit calibration spectra were taken on the first 
and third night immediately after each of the three object spectra. On the 



^ Based on observations made with ESO telescopes at the Par anal Observatory 
under program IDs 68.A-0106 and 70.A-0017 
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second night there is no long-sht cahbration spectrum after the second object 
spectrum so only two long -sht spectra ex i st for that night. This appears to 
contradict the statement in llvanchik et al.l (120051 ) which states that calibration 
spectra occurred before and after each object spectrum but it is possible that 
there are other calibration images that the archive did not associate with the 
QSO integrations. The Q0347-383 observations were taken at a grating setting 
that put the center wavelength at 4300 A and the Q0405-443 observations at 
a center wavelength of 3900 A. We will refer to these as the 430 and 390 grat- 
ing settings in the following. We use only the blue channel spectra since they 
contain the spectral area of interest for the measurement of fi. Table [T] gives 
some of the observing parameters for the long-slit calibration lamp images. 



2.1 Suggested Observing Procedure for Cahbration Spectra 



For very accurate wavelength calibration we suggest an observing sequence 
that takes several (more than five) long-slit calibration line spectra immedi- 
ately before and immediately after each spectral image of the observed object. 
The primary purpose of the multiple calibration line images is to produce a 
meaningful median image that is free of cosmic ray events rather than an 
increase in the signal to noise ratio of the image. The time pacing of the cal- 
ibration images also tracks drifts in the wavelength calibration due to any 
contributing factor. The small amount of extra time devoted to this observing 
sequence can make the production of cosmic ray free images straight forward. 
Observing the calibration lamp immediately before and immediately after an 
object observation provides an accurate monitor of any time variation in the 
wavelength calibration. If the object observations are spaced close together 
one set of calibration lamp exposures can serve as both the before and after 
function for intermediate object exposures. 



3 The Pipeline Data reduction 



All files associated with the observations of the two QSOs were downloaded 
from the ESO data archive. The MIDAS based UVES pipeline reduction soft- 
ware was also downloaded and installed in a LINUX based operating system. 
All of the data were processed with the pipeline software which removes the 
bias and dark counts. The bias and dark corrected files are designated by the 
the suffix b_.bdf in their file names. All of the image files are converted to FITS 
format files using the MIDAS OUTDISK/FITS procedure. Our analysis starts 
with the two dimensional FITS format images of the long-slit calibration lamp 
images and long-slit fiat-field images that have been bias and dark subtracted 
but not fiat fielded. These files are identified by the terms LAMP, WAVE and 
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LAMP,FLAT in the OBJECT parameter in the file header. We also utilize 
the short slit lamp flats identified by LAMP,ORDERDEF for identifying the 
location of the orders. A pipeline produced table of x and y line positions is 
used to provide a first guess and the wavelength solution. In this and in all 
other references to x and y, the x direction is along the dispersion direction 
and the y direction is perpendicular to the dispersion direction as seen in the 
pipeline rotated images. The order traces are generally along the x direction 
but are curved and tilted upward. The pipeline produced intermediate product 
spectra are used to produce the final spectrum. We do not use the primary 
output of the pipeline which are spectra interpolated to an equally spaced 
wavelength grid. 



4 Wavelength Calibration 

All procedures are written in IDlH] code. In the following we will write IDL 
provided procedure names in capital italic letters and procedures written by 
the authors using IDL code in lower case italics. After further refinement 
the authors intend to make the reduction code available to the public, hence 
the use of the procedure names in the flow charts. Following the production 
of the pipeline products the wavelength calibration proceeds in two semi- 
independent phases. The first phase, called the first pass, uses the long-slit 
bias corrected images produced by the pipeline reduction. In the first pass 
the long-slit calibration images for a given setting are median combined to 
form a "master" long-slit calibration image that is free of cosmic ray hits. 
A first guess at the wavelength solution based on the line positions output 
by the pipeline provides the initial wavelength solution for each order. The 
solution is then iterated interactively and a set of "good" calibration lines is 
identified. This step involves some subjective judgment on the definition of 
good that is discussed below. In the second phase, called the second pass, 
the exact spatial locations of the good lines are determined. Each long-slit 
calibration image is shifted to a common position and a new median image 
is produced. The shifts are typically on the order of a few hundredths of a 
pixel. The wavelength solution is again interactively iterated to provide the 
final solution for each order. The primary output is a high signal-to-noise ratio 
master long-slit calibration lamp image and a vacuum wavelength solution for 
each order. The method assumes that the differences between two calibration 
lamp images at the same grating setting are small in order to form a useful 
median image at the beginning of the first pass data reduction. Once the 
master long-slit calibration line image has been produced, however, the actual 
shifts of subsequent images can be accurately determined from the master. 



^ IDL stands for Interactive Data Language registered by Research Systems Inc. 
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An overview flow chart of this process is shown in Figure [H The individual 
procedures of the analysis are described in subsequent sections. 

First pass wavelength solution 

The flrst phase of the method, after producing the pipeline products, is to make 
a flrst pass at the wavelength solution which will be improved in the second 
pass. This step can be eliminated in subsequent measurements at the same 
grating setting once a master long-slit calibration line image and wavelength 
solution have been established. A flow diagram of the flrst pass wavelength 
calibration procedures is shown in Figure [21 The names in the boxes that end 
in .pro are simply the names of the IDL based procedures that were developed 
for this process. 

The flrst pass starts by sorting all of the two dimensional images to determine 
their type by the entry into the OBJECT parameter of the image headers. All 
of the long slit calibration lamp images are median combined to eliminate the 
cosmic ray events. There is no attempt in this flrst pass to adjust the images 
for small shifts in the location of the lines. The positions of the lines typically 
had shifts of a few hundredths of a pixel over a three-night set of images. 
The 390 grating setting median contains 6 images and the 430 setting image 
contains 8 images. Visual inspection of the images did not detect any residual 
cosmic ray hits. 

4-. 1.1 Order definition 

The location of each of the spectral orders must be determined to start the 
spectrum fltting. The short slit flat lamp images are used to determine the 
order positions. The flrst step is to median the images together to elliminate 
cosmic ray events. One order deflnition image was taken in each night there- 
fore the medians contain three images which is a minimum, but certainly not 
optimal, number of images to form a median. Figure [3] shows the 1500 x 1024 
binned pixel median order-deflnition-image for the 390 grating setting. Fig- 
ure [3] shows that the intensity of the order deflnition image varies strongly over 
the image in both the horizontal and vertical directions. This is a challenge to 
an automated order flnding algorithm. The intensity variation of the spectra 
is similar and results in the flrst 100 pixels of the spectrum being rejected due 
to a poor signal-to-noise ratio. The orders are curved and are slanted upward 
going from left to right. 

Next a flrst solution to the position of the orders is found. At this point the 
procedure becomes partially UVES speciflc. The procedure knows that the 
pipeline has rotated the images so that the orders run horizontally from left 
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to right with a shght upward tih. It therefore looks for new orders appearing 
from the bottom of the image and knows that it should look for orders along 
the vertical direction. An efficiency map is formed by smoothing the order 
image with a 101 by 101 box four times. The search for orders starts 100 pixels 
past the left edge of the order image. The region past 100 pixels is binned in 
groups of 10 columns each and the groups are summed to form single columns 
of higher signal to noise for both the original order image and the smoothed 
image. The derivative is taken of each order column and the order position 
is defined by the S-shaped curve produced when an order is encountered. 
The order position is where the value of the derivative changes sign between 
positive and negative. The derivative must have at least two positive values 
followed by two negative values to be considered an order position. The flux 
at that position in the summed order column is compared to the flux in the 
smoothed column to check against noise that could produce the same pattern 
of two positive values followed by two negative values in the derivative. If an 
order is found below the lowest existing order a new order is declared and it 
becomes the lowest existing order. 

After all of the order points have been found the orders are fit in x and y 
with a quadratic fit. The fit is displayed on the computer screen overlaying 
the order image to check for any bad fits. The output is an array which has 
the quadratic coefficients plus the pixel numbers of the starting and ending 
pixels where the order is valid. Orders that run off the top of the image have 
ending pixel values less than the 1500 pixel extent of the image. This array is 
used as the starting point of the final step in the order definition. 

After the initial order definition step the y position of each order for each 
pixel column is defined. At each order position defined by the quadratic fit 
a Gaussian fit is found in an 11 pixel column centered on the pixel position 
given by the fit. The IDL procedure GA USSFIT fits a Gaussian plus a constant 
offset to the column to determine the y position of the peak. After the x and 
y position for each order is determined by Gaussian fitting each order is fit 
with a six term Legendre polynomial fit using the IDL procedure SVDFIT 
with the Legendre option set. The Legendre polynomial order fit gives the y 
position in pixels of each order for every pixel column. This is then taken as 
the vertical center position of the order at each horizontal pixel position in 
the long slit calibration lamp images in the following analysis. 



4.1.2 Flat field 

Although not absolutely critical to the wavelength calibration, the next step in 
the analysis is the production of the flat-fleld. This is identical to the flat fleld 
used in the object spectrum analysis, therefore it is also used in extracting the 
calibration lamp spectrum for consistency. The fiat-field is constructed from 
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the median of all of the long-slit flat-field images taken at the desired grating 
setting. The procedure first determines the length of the slit by looking at 
the center of the middle order. It determines along a vertical line where the 
intensity of the image falls to 1/2 of the maximum and defines the height 
between those points as the length of the slit. The procedure produces a fiat 
correction image which is the inverse of the long-slit fiat lamp image in the 
region of each order defined by the length of the slit and the beginning and 
ending of the order as defined in § 14.1.11 The image is divided by the median 
of all of the non-zero pixels in the fiat correction image to make the median of 
the non-zero pixels equal to one. Correction of the object spectra is achieved 
by multiplying the object spectrum image by the fiat correction image. 

4- 1-3 Stray light 

The long-slit calibration lamp images have low level background stray light 
from nonspecular reflections, grating imperfections and other sources. The 
next step in the analysis removes the smooth component of the stray back- 
ground light to first order. The photon noise from the background light is, 
of course, still in the stray-light-corrected image. To start the procedure the 
median long-slit calibration lamp image is multiplied by a mask that is zero in 
the regions that contain the calibration lamp spectra and 1 in the inter-order 
regions. We assume that the stray light observed in the inter-order region ex- 
tends smoothly into the regions containing the orders. The masked image is 
divided into regions of 10 columns each. Each region of 10 columns is median 
combined in the horizontal direction to form single columns of higher signal-to- 
noise ratio. A third order polynomial is fit along the vertical direction of each 
of the columns using only the non-zero inter-order regions of the columns. This 
produces a set of measures of the stray light in the vertical direction spaced 
every 10 columns in the middle of each 10 column region. This results in a 
stray light image that has more resolution in the horizontal direction than in 
the vertical. This was done because the variation of the stray light appears 
to be stronger in the horizontal direction than in the vertical. The final stray 
light image is formed by interpolating the solution in the horizontal direction 
for each row of the image. This smooth distribution is then subtracted from 
the median long-slit calibration lamp image. 

4-1-4 Extraction of the spectrum 

The spectrum extraction procedure produces a two dimensional image of each 
order. The height of the image is the slit length in pixels and its length is the 
length of the spectral image in pixels. The total output is a four dimensional 
array with the dimensions spec [nx,norders, 6, length] where nx is the number 
of pixels in the x direction, norders is the number of usable orders and length 
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is the length of the sht in pixels. For the 6 inputs in the third index the first is 
the wavelength, which at this stage is just the pixel number, the second is the 
flux, the third the noise, and the other 3 are populated later in other spectral 
analysis procedures. The first step in the procedure is to multiply the total 
spectral image by the fiat-field correction image described in § I4.1.2[ 

The procedure next extracts the two dimension spectral image of each order. 
The center pixel in the vertical direction for each order is determined by the 
order solution discussed in § I4.1.1[ The IDL procedure ROUND finds the 
pixel number and the pixels above and below the center pixel are extracted 
according the the value of the slit length in pixels. The flux is in ADUs (Analog 
to Digital converter Units). The spectral image of the fifth order in the array 
is spec[0:nx-l,4,l,0:length-l] in the IDL array. In IDL the nomenclature a:b 
means all entries with indexes a through b. The first index in IDL is rather 
than 1 which accounts for the length- 1 in the preceding and why the fifth order 
has index 4. The noise for each pixel is computed as the root mean square of the 
photon and read noise. The flux in ADUs per second is converted to electrons 
using the conversion factor in the header and the fiat-field correction. The read 
noise is also listed in the header. The calculated noise in electrons is the square 
root of the sum of the total number of electrons plus the square of the read 
noise. The noise in electrons is converted back into ADUs and re-corrected for 
the fiat-field. The noise image is in the same format as the flux except that 
the third index in the extraction array is 2 rather than 1. 



4- 1-5 Spectrum assembly 

The next step is the creation of a one dimensional spectrum for each order 
from the individual strips of spectra created by the extraction. Inspection 
of the long-slit calibration line image indicates that the intensity along the 
length of the slit is quite uniform and that the image of the slit is vertical 
with no apparent deviation over the entire field of view. For these reasons the 
assembled one dimensional spectrum is simply the sum of all of the pixels in 
the vertical direction. If there is any fractional pixel tilt in any part of the 
field of view it would result in a slightly lower resolution but the center point 
of the line would still be the same. A comparison of the object spectra with 
the order definition positions confirmed that the object spectra maxima lie 
at the middle of the slit length. The noise is simply the square root of the 
sum of the squares of the noise determined for each pixel during the spectrum 
extraction. After this stage we have individual one dimensional spectra for 
each order. The physical orders are 141 to 102 for the 390 grating setting 
and 127 to 95 for the 430 grating setting, however only a subset of the orders 
had object spectra appropriate for the determination of the value of /z at 
the redshifts of the absorption line systems. We list the order numbers in 
inverse order since the image orientation returned by the pipeline has the 
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order number decreasing in the vertical direction. The spectrum data cube 
now has the dimensions spec[nx,nord,6] since the previous array was summed 
along its fourth dimension. Note that nord equal to zero corresponds to order 
141 and 127 for the 390 and 430 grating settings. A value of nord equal to 1 
corresponds to orders 140 and 126 respectively etc. 



4.1.6 Preliminary wavelength assignment 

The UVES pipeline reduction produces a table of the x and y positions for 
calibration lines during the pipeline reduction along with the air wavelengths 
of the line and the order number. The table is labeled IxxxBLUE.tbl for the 
blue channel where xxx is the grating setting such as 390. We use the pro- 
cedure lineguess.pro on the output of this table to produce a preliminary 
wavelength solution for each order. The first step is th e conversion of the air 



wavele ngths to vacuum wavelengths by the formu l as of lEdlenl (119661 ) as given 



in the ICRC Handbook of Chemistry and Physics! (119941 ). The lineguess pro- 



cedure uses the IDL procedure SVDFIT.PRO in double precision with the 
Legendre polynomial option set to produce a six term Legendre polynomial 
fit to the wavelengths. The outputs of the procedure are six coefficients of 
the first six Legendre polynomials for each order. These coefficients are the 
inputs to the procedure waveputlml .pro which converts the pixel numbers 
in the wavelength array of the spectrum to wavelengths. The wavelength for 
each of the nx pixels for an order nord in the output data cube is given by 
spec[0:nx-l,nord-l,0]. See § 14.1.101 for more details of the wavelength fitting. 



4.1.7 Thorium and Argon line lists 



An important part of the wavelength calibration improvement over the pipeline 
calibration is the use of newly pub lished line list s . The Tho rium and Ar- 



gon li n e lists came from tw' o sources, iLovis fc Pepd (120071 ) and iMurphy et al. 



(120071 ). iLovis fc Pepd (120071 ) list all of the Thorium- Argon lines between 3785.6 
and 6914.4 A found in a dedicated observati on of Thoriurn - Argon lines with 
the HARP spectrometer. The line list from iMurphy et al.l (120071 ) lists lines 
between 3033.5 and 1 0506.5 A. This list o nly includes lines that have been 
deemed acceptable in IMurphy et al.l (120071 ) for accurate wavelength calibra- 
tion. Our lists were constructed from files kindly sent by Michael Murphy prior 
to publication, however, the files are now available from websit es listed in the 
public ations. Our final list of li nes contains al l of the lines in the lLovis fc Pepe 
(I2OO7I ) list plus the lines from iMurphv et al.l (120071 ) that lie at shorter wave- 
lengths than the Lovis fc Pepe ( 20071 ) list. The data from the lists for each line 
are entered into an IDL data structure. An IDL data structure is a list of pa- 
rameters that can be a mixed combination of numbers, text strings or arrays, 
each with a different format. The structures for all of the identified lines are 
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identical in format and the entire set of lines forms an IDL structure array for 
easy access by the IDL based analysis procedures. The data items for each line 
in the IDL structure are: the wavenumber in cm~^; the vacuum wavelength 
in A, the standard pressure and temper ature (STP) a ir wave length, the pre- 
viously determined w avelength listed by lLovis fc Pepd (120071 ). the error from 



Lovis fc Pepd (120071 ). the intensity, the identificati on, and a "Murphy" flag. 



The "Murphy" flag is the letter M if it is listed by iMurphv et all (1200711 as a 
good line and the letter N if it is not. Note that the iLovis fc Pepd (120071) list 



contains all of the observed lines in its wavelength range while iMurphy et al 



( 120(371) only list the goo d lines therefore all of the lines outside the range of 
the iLovis fc Pepd (120071 ) lists are "Murphy" lines. 



4-. 1.8 First determination of the line positions and identifications 

The next step in the first pass through the wavelength calibration is the lo- 
cation of the calibration line positions in the median long-slit calibration line 
image and a preliminary identification of the proper Thorium or Argon line 
for the line position from the line list described in § 14.1.71 

The line identification procedure for a particular order first determines the 
median flux in the order which is assumed to be the continuum emission 
from the calibration lamp and subtracts that from the order flux. Possible line 
positions are identified by taking the derivative of the flux versus pixel position 
with the IDL procedure DERIV and looking for the positive to negative zero 
crossing that the positive and negative slopes of an emission line create in the 
derivative. A legitimate zero crossing requires two positive derivative values 
followed by two negative derivative values. The procedure does several checks 
on the legitimacy of the position. The flux at the position must be greater 
than a preset threshold which is equal to the median previously subtracted 
from order flux. If the line flux is only positive for one pixel or if the line 
appears on the side of another line it is rejected. 

After identification via the derivative method a Gaussian fit is done at the 
positions that survived the checks described above to find the line position with 
fractional pixel accuracy. This is a straight Gaussian fit with no continuum as 
the continuum was subtracted earlier. This is done with a slightly modified 
version of the IDL procedure GAUSSFIT called gaussfitstat that returns the 
status flag of the fitting procedure as well as the parameters of the Gaussian 
fit. Any position with a bad status flag, indicating an internal IDL error such 
as underflow or overflow, is rejected. The position found by the Gaussian flt is 
placed in a new line structure along with the coefficients of the Gaussian fit. 
Unlike the Th/Ar line list the new structure array only contains structures for 
the lines found in the median long-slit calibration line image. The structure 
for each line has 22 different entries which eventually include the information 
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from the Th/Ar line lists, the "Murphy" flag, the spectral order and other 
information for record keeping. One of the entries is a "use" flag on whether 
it should be used in the flnal wavelength solution. Not all lines with positive 
"Murphy" flags were used in the solution. 

After all of the line positions have been determined a preliminary assignment of 
line identiflcations is made. The preliminary wavelengths are matched against 
an atomic line list described above in § I4.1.7[ The line in the list with the 
closest wavelength to the measured wavelength of the line is assigned as the 
preliminary identiflcation. There are 3751 identifled lines for the 390 grating 
setting and 3087 hues for the 430 setting. Most of these lines are not be used 
in the wavelength calibration 



4- 1-9 Using vacuum wavelengths 

The wavelength calibration is performed using vacuum wavelengths. This 
means that the wavelength solution produced assigns the proper vacuum wave- 
length for each Th/Ar line as the appropriate wavelength for that pixel posi- 
tion in each order. The observations, however, were in conditions nearer to STP 
than vacuum. The wavelength calibration solution simply gives the vacuum 
wavelength as a function of pixel position. As such it is only valid for observa- 
tions carried out under similar atmospheric conditions to those present when 
the calibration observations were taken. This would be true whether we used 
STP wavelengths, current atmosphere wavelengths or vacuum wavelengths. 



The iLovis fc Pepd (120071 ) measurements were made in vacuum therefore we 
have chosen vacuum wavelengths as our standard. The shifts in temperature 
for any night of observation were 0.1 degrees Centigrade or less and the pres- 
sure shifts were 0.7 mm. Hg or less during a night as shown in Table [U At 
4000 A the rate of change of air wavelength is less than 4 x 10"'^ A per degree 
Centigrade and slightly more than 1 x 10~^ A per mm. Hg. These are changes 
on the order of 1 part in 10^ or less for all nights and are 1/3 or less of the 
accuracy of our wavelength solution. Since shifts are calculated on a night to 
night basis and on an order by order basis, any shifts between nights 
counted for in the zero point calculation for the wavelengths of each order for 
a night. Future improvements, such as higher resolution spectra, will require 
a temperature and pressure correction in the procedure. 



4.1.10 Refining the wavelength solutions 

The flnal step in the flrst pass through the wavelength calibration is a re- 
flnement of the wavelength solution. This is accomplished with an interactive 
IDL procedure called waverefine^uves in Figure [H The procedure provides a 
method of improving the preliminary 6 term Legendre polynomial wavelength 
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solution. Wavelength solutions employing different polynomials and Legen- 
dre polynomial fits of orders 4 through 7 were tried initially. The Legendre 
polynomial fits appeared to be the most stable, without high deviations from 
acceptable fits at the end points of the solutions. The quality of fits improved 
up to the 6 order solution used in this work. The addition of an additional term 
to make a 7 order fit did not improve the fits and in some cases made them 
worse. This choice is probably unique to the UVES blue channel and may be 
different for other spectrometers with different sources of image distortion. 



For each order the procedure initially deems all identified lines as not usable. 
As each section of the spectrum for an order is displayed on the monitor the 
user indicates which lines should be included in the wavelength solution. A 
vertical line for each identified line is displayed at the wavelength listed in the 
line list with a height that is proportional to the strength in the line list. This 
helps the user match the observed spectrum to the appropria te lines in the 



l ist. If the line is one of the vetted acceptable lines identified in iMurphy et al 



( 120071 ) then a letter M is printed above the line position. During this first 
pass the majority of the lines are rejected based on a poor signal-to-noise 
ratio, asymmetry due to blending or obvious double lines. Each order contains 
between 80 and 100 lines from the list described above. Usually about 25 lines 
are selected as usable. The choices of usable lines are subjective but most of 
the rejectio ns are fairly obvious. This rejection ratio is similar to the ratio 
described in iMurphv et al.l (120071 ). Our selection of lines is independent of the 



Murphy et al.l (120071 ) selections in the following manner. Our inital selection 



and rejection of lines was done without reference to which lines were Murphy 
lines. This led to the result that some lines that were not Murphy lines are 
incl uded and seve r al Mu rphy lines were rejected. For lines that are bluer than 
the iLovis fc Pepd (120071 ) list only the Murphy lines were available but many 
of those were rejected. The statistics for the 390 grating setting are that there 
are 2209 Murphy lines in the spectral region of the observations. Of these 
we chose 518 to use in our analysis plus 2 that are not Murphy lines. For 
the 430 grating setting there are 1462 Murphy lines in the spectral region of 
which we chose 527 plus again 2 lines that are not Murphy lines. For obvious 
reasons most of the selected lines are Murphy lines due to similar selection 
criteria. Almost all of the rejections of Murphy lines in our list are due to a 
more conservative cut on line crowding. This is a key difference in using the 
2x2 binned pixels rather than the individual CCD pixels. Another difference 
is that our selection is order by order for each of the grating settings. A line 
that is selected in one order may be rejected in an adjacent order if it falls in 
a poor signal to noise region. 



After the first pass through the order the procedure produces the Legendre 
polynomial fit to the lines marked usable and displays the deviation of the 
measured lines from the fit. Typically there will be a few extreme outliers. 
The program allows the user to add or remove lines from the usable list one at 
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a time. The new solution is displayed after each addition or removal. This is 
the most subjective part of the procedure and the most scientifically dubious. 
The rejection of outlying data points is risky. The assumption here is that the 
true wavelength solution should be smoothly varying without any discontinu- 
ities or high-frequency variation. Practice with the procedure leads to better 
identification of the outliers that are causing the errors in the solution. Often 
it is not the furthest outlier that is causing a problem. Removal of a single out- 
lier that produces discontinuous break between positive and negative outliers 
often greatly reduces the average error. 

Another subjective aspect to using the procedure is that the final set of de- 
viations from the solution should be "random". A solution that produces a 
sharp break in the sign of the outliers or that has a systematic trend in the 
location of the outliers is probably wrong. After the solutions are adjusted 
typically 20 of the original 25 hues are left as usable. After all of the orders 
have their solutions the procedure is repeated several times to improve the fit. 
This allows some lines that may have been erroneously rejected in a previous 
pass to be returned to the usable line list and may result in rejection of other 
lines. 



4.1.11 First pass data products 

The first pass data products are: the first pass wavelength solutions which is a 
set of the 6 coefficients of the first 6 Legendre polynomials for each order that 
represent the wavelength solution for that order in terms of A versus pixel; 
a median long-slit calibration image, the first pass one dimensional order by 
order spectrum; and the long slit calibration line structure array with the use 
fiag set positive for all of the lines that were deemed usable in the first pass. 
All of these products are updated in the second pass through the wavelength 
solution. The flat correction image, order deflnition solutions and the length 
of the slit are first pass data products that are not updated in the second pass. 



4.2 Second pass wavelength solution 

The second pass to the wavelength calibration recognizes that the individual 
long-slit calibration lamp images taken over the course of the observations 
are not perfectly ahgned. It computes the offsets between the images using 
the lines that were declared usable in the first pass. The procedure shifts the 
images to a common registration and produces a new median image which 
becomes the master long-slit calibration lamp image. It then follows most of 
the steps used in the first pass to produce a new wavelength solution for the 
master long-slit calibration lamp image. It also preserves the shifts calculated 



15 



for the individual images to provide a baseline to compute the shifts appropri- 
ate for the source spectrum images. A flow chart of the second pass procedures 
is given in Figure |H 



4.2. 1 Calculating the long-slit calibration image ojfsets 

The procedure starts by reading in the list of long-slit calibration lamp im- 
ages and the line structure array generated in the first pass. For each line 
marked usable the procedure make Jong caLmask finds the corner coordinates 
in pixels of a 7 by 21 rectangle in x and y centered on the line position. These 
small sub-images are the regions that will be cross correlated to find the shifts 
between the images. The first long-slit calibration lamp image is arbitrarily 
chosen as the reference image. For each line in each image the small two di- 
mensional image is summed in the y cross dispersion direction to form a 7 pixel 
one dimensional spectrum. The one dimensional spectra are interpolated to a 
grid that has a 0.01 pixel spacing and cross correlated with the corresponding 
interpolated spectrum in the reference spectrum. The lag producing the max- 
imum cross correlation signal is selected as the shift for that line relative to 
the reference image. After all of the shifts have been calculated the shift for 
each long-slit calibration line image is taken as the median of the shifts for all 
of the usable lines in the image. Histograms of the shifts for each image are 
also plotted to check for anomalous distributions of shifts. The shifts for the 
390 and 430 images are given is Table [H 



Once the shifts have been calculated each image is shifted to the position of 
the zero shift image by interpolation. The IDL procedure INTERPOLATE 
is used with the value of the CUBIC parameter set to -0.5 and the keyword 
GRID set to true. The value of CUBIC is the cubic polynomial interpolation 
parameter. It is the value of a in the following cubic polynomials r(x) which 
are used t o calculate the contribution t o the interpolated value from adjacent 



pixels, see iPark &: Schowengerdtl (119831 ) . 



r{x) = (« + 2) I X 1^ -{a + 3) | x ^ +1 
r(x) = a I X P —5a I x P +80; I x I —4a 



I X |< 1 
Kl X |< 2 



r(x) =0 otherwise 



The shifted images are median combined to produce a cosmic ray free long-slit 
calibration lamp image. The stray light removal procedure described earlier in 
§ I4.1.3l is performed on the median image and the result is stored as the master 
long-slit calibration lamp image. Future long-slit calibration lamp images can 
be added to the master image by repeating the processing with the full set of 
images. 
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The second pass spectrum extraction and spectrum assembly are identical to 
the first pass processing. The wavelength solution found from the first pass is 
put into the spectrum. The IDL procedure firstlineposition.pro again assembles 
the detected lines into a line structure array with the new parameters from 
the master long calibration lamp spectrum. Rather than starting from scratch 
the line usage flags from the first pass are transferred from the first pass line 
structure array to the second pass array by the procedure lineinfoexchange.pro. 

The wavelength solution is again refined in iterative passes as described in 
§ I4.1.10[ The result of this procedure is the master wavelength solution. Out 
of the 3751 identified lines in the 390 grating setting 686 lines were used in 
the wavelength solution. Only 4 are not in the vetted Murphy list. 577 lines 
out of the 3087 identified hues were used in the 430 grating setting wavelength 
solution, 3 of which were not in the Murphy list. The difference between these 
numbers and the numbers in § 14.1.101 is that here a line is counted more than 
once if it appears in more than one order. The use rate of the identified lines 
iss 18% for the 390 setti ng and 19% for the 43 setting. This is lower than the 



retention fraction in the iMurphy et al.l (120071 ) list. 



4-2.2 Final wavelength calibration data products 

The final data products for each grating setting are stored for use in determin- 
ing the wavelength solution for the object spectra. They include the master 
long-slit calibration lamp image, the master one dimensional order by order 
spectrum, the master wavelength solution for that spectrum and the associ- 
ated structure array of line parameters. Also stored are the wavelength shifts 
in pixels for each of the individual long-slit calibration lamp images and the 
times that they were taken. The last data sets are the median fiat correction 
image, the order definition solutions and the length of the long-slit in pixels 
from the first pass solutions. A list of the UVES Thorium and Argon lines 
used for the 390 and 430 grating settings wavelength solutions is included in 
the electronic version of this paper. These list s are particular to th e grating 



settings as opposed to the vetted line lists in iMurphy et al.l (120071 ). This is 
because a line that has a good signal-to-noise ratio in one order may fall in a 
poor signal-to- noise ratio area in an adjacent order. 



4-3 The form of the wavelength solutions 



The wavelength solutions are given by 

5 
i=0 
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where Li are the Legendre polynomials and the Cj are the coefficients for each 
order that are the products of refinewave-uves.pro. The Legendre polynomials 
are only valid in the range between -1 and 1 therefore the variable x is related 
to the pixel number p by 

P- 750. 

X = 2 

750. ^ ^ 



where p ranges between and 1499 for the 1500 pixel width of the blue channel. 
This means that the largest value of x is 0.99866666 rather than 1 but that is of 
no consequence since each value of x is unique, equally spaced and within the 
valid range of the Legendre polynomials. Given the nature of the Legendre 
polynomials the coefficient Co is the wavelength in A of the middle of the 
order and Ci is the linear dispersion in A per value of x. The higher order 
coefficients describe the deviation from linearity of the wavelength solution. 
Tables [2] and [3] give the coefficients found for the master long-slit spectra 
at the two grating settings. Since all of the long-slit calibration images were 
shifted to the position of the ffist calibration image for each run the solutions 
presented here represent the best wavelength solution for the ffist calibration 
image. Since the shifts are small between images these solutions can be used 
as starting points for other spectra at those grating settings. 



4.4 Accuracy of the solutions 



We can estimate the accuracy of the wavelength solution by looking at how 
accurately the 6 term Legendre polynomial solution reproduces the measured 
vacuum wavelengths of the calibration lines used in producing the solution. 
The statistical look at the fits is appropriate since for most orders the number 
of lines used in the solution is 3 to 4 times the number polynomial terms. 
The statistics of the fits for all of the lines gives a standard deviation on the 
fractional error AA/A of 2.6 X 10""^ for the 390 grating setting and 2.2 x lO"'^ 
for the 430 grating setting. A lis ting of the fracti o nal e rrors for each order 
used in the program described by [Thompson et al.l (120081 ) is given in table HI 



The distribution of errors is shown in Figures [5] and [6] where the fractional 
differences between the calculated and listed wavelengths are shown in his- 
tograms with a bin size of 1 x 10"''. The histograms are reasonably symmetric 
about zero and the la values corres pond to roughly 75 m s~^. Comparison 
with Figure 8 in lMurphy et al.l (120071 ) indicates that the wavelength errors are 
a f actor of 2 less t han t he pipeline values but somewhat more than achieved 
by iMurphy et al.l (120071 ) for their data sets which contain more calibration 
line images than this work, a narrower slit width (0.6 rather than 0.8 arc 
seconds) and unbinned images. In both cases the line selection is the primary 
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component in improving the accuracy of the cahbration. The data sets we use, 
however, are the cahbration hne spectra taken at the time of the observations 
we want to cahbrate and better represent the wavelength cahbration for those 
objects. 



One of the primary points of iMurphy et aLl (120071 ) is that the pipehne UVES 
wavelength calibration has significant systematic residuals as shown in their 
Figure 9 and that this residual can give false indications of variation in the 
measurement of the fundamental constant a. The same systematic e rror resid- 



ual ca n affect the measurement of as well but in a different manner (IThompson et al. 



20081 ) than for a. Since the primary purpose of this paper is to document the 
wavelength calibration for the determination of /i at e arly times in the uni- 
verse we have carried out an analysis similar to that of iMurphy et al.l (120071 ) 
for our wavelength calibration. Figures [7] and [8] show the scatter of all of the 
calibration line fits that are used in the determination of [i at both the 390 
an d 430 grating s etting s. This figure can be directly compared to Figure 9 



m 



Murphy et al.l (120071 ). The solid black line is the mean residual in 100 A 



bins. It is clear that the large scale residuals present in the UVES pipeline 
calibration are not present in the wavelength calibration presented here. The 
minimum near 3550 A in figure [7] is not well understood but may be simply a 
statistical anomaly. 



5 Addition and wavelength solution of the individual object spec- 
tra 



Once the wavelength calibration has been performed it can be applied to the 
individual spectra of the objects produced by the pipeline. A typical script 
for the pipeline reduction is given in Appendix |A] The standard output of the 
UVES pipeline are spectra that have been interpolated onto a constant incre- 
ment wavelength grid. All orders have been combined to produce a continuous 
spectrum with no order breaks. The details of this proc ess are not im medi- 



ately obvious from the pipeline documentat ion dB alles ter et al.l. 120041) . The 



results of the optimal extraction described in lBallester et al.l (120041 ) are, how- 
ever, available as an intermediate product. These files for the blue channel are 
labeled fxb_UVES*_b.fits, where the * represents a date and time tag along 
with a 3 digit identifier. These are order by order one dimensional optimal 
extractions of the spectra. 
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5.1 The wavelength solutions for the object spectra 



The object spectra wavelength solutions start with a determination of the 
shift in wavelength of the object spectrum relative to the master wavelength 
solution. This is done by fitting the shifts of the individual long-slit calibration 
lamp images versus time for each night of observing. For nights with just two 
images this is a linear fit. Nights with 3 images are fit with a quadratic fit. 
The shifts for the object spectra are determined from the shift versus time fit 
function by the midpoint in time of the integration. The shifts are similar to 
the shifts for the calibration spectra shown in Table [T] which are on the order 
0.01-0.02 pixels. The master solution is interpolated to the shifted positions 
found above and placed into the wavelength section of the spectrum array. This 
wavelength solution is the appropriate vacuum wavelengths for a source at rest 
relative to the spectrometer and therefore has to be corrected for barycentric 
velocity. The appropriate barycentric velocity is calculated for the midpoint of 
the observation and the wavelength solution is then corrected to account for 
the barycentric velocity of the earth at the location of the VLT. At this point 
in the analysis the center points of the pixels in the different object spectra 
have different wavelength solutions. 

In preparation for addition of the spectra the spectra are all shifted to the 
wavelength solution of the master wavelength solution by interpolation. The 
interpolation shift is determined from the shift relative to the master solu- 
tion calculated above. This is not done initially since the various spectra have 
slightly different barycentric velocity contributions. As with the long-slit cal- 
ibration spectra that bracket the observations, the shifts from the master 
solution is generally on the order of a few hundredths of a pixel. The spectra 
now all have the same pixel wavelength solution and can be added together. 



5.2 Addition of the spectra 



Since the spectra now all have the same wavelength per pixel solution they 
can be directly added together for each order. The 9 spectra for each object 
are combined to produce both the median and the mean for each pixel. The 
median better rejects any residual cosmic ray hits in the spe ctra and provides 



the pr imary spectrum for the analysis of the value of in [Thompson et al 



(120081 ) . The mean spectrum is also analyzed to provide a check on the median 
spectrum. Other "optimal" additions methods were utilized that weight the 
pixel flux by its signal to noise ratio. These methods did not produce quantifi- 
bly better results than the median spectra. This is most probably due to the 
remarkable uniformity of the spectra between observations of the same ob- 
ject. For this reason we utilize the median spectra in our analysis. Data sets 
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with significantly varying signal to noise ratios between spectra will certainly 
benefit from more optimal spectrum addition methods. 



6 Conclusions 



The measurement of fundamental constants through astronomical spectra re- 
quires the best possible wavelength calibration that can be achieved. This re- 
quires both observing and data reduction procedures beyond what is practical 
for most standard observation and data reduction procedures. The procedures 
described in this paper have achieved a la fractional AA/A wavelength accu- 
racy of better than 3 x 10"'' for UVES spectra which is a factor of 2 better 
than the accuracy of the UVES pipeline reduction. More impor tantly there is 



no det ectable systematic residual wander as was discovered by iMurphy et al 



( 120071 ) in the pipeline wavelength calibration. 



The reduction methods are described in detail so that they can be repeated 
by other researchers and a suggested observing sequence was presented that 
can improve the accuracy of both the wavelength calibration and its transfer 
to the object spectra. This paper serves as a detailed description of the cali- 
bration procedures used in the determination of the fundamental constant ^ 
through the molecular hydrogen absorption lines in the spectra of two QSOs 



Thompson et al.l ((20081). 
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Grating Setting 



Grating Setting 





390 








430 






MJD Date 


shift 


Temp. 


Press 


MJD Date 


shift 


Temp. 


Press 


JD-2.4E6 


pixels 


C° 


mm. Hg 


JD-2.4E6 


pixels 


C 


mm. Hg 


52643.08644 


0.00 


14.8 


743.6 


52282.08493 


0.00 


11.8 


745.6 


52643.14617 


-0.06 


14.8 


743.6 


52282.13882 


-0.02 


11.8 


745.7 


52643.20522 


-0.05 


14.7 


743.0 


- 


- 


- 


- 


52644.09176 


-0.02 


14.5 


743.3 


52283.08329 


0.03 


11.2 


746.4 


52644.21381 


-0.03 


14.4 


742.2 


52283.13777 


0.04 


11.2 


746.4 


52645.09255 


0.00 


13.5 


744.2 


52284.08694 


0.01 


11.5 


745.6 


52645.15520 


0.01 


13.5 


743.9 


52284.14110 


0.02 


11.5 


745.7 


52645.21402 


-0.02 


13.5 


743.2 











Table 1 

Dates, shifts in position in pixels, temperature and pressure for the 390 and 430 
grating setting long-slit calibration lamp images 
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order 




Co 




Ci 




C2 




Cs 


C4 


Cs 


139 


3, 


.3578451E+03 


2, 


.6241278E+01 


-2 


.2600063E+00 


4, 


.0406710E-02 


-1.5302832E-02 


8.2013765E-03 


138 


3, 


.3821772E+03 


2, 


.6446218E+01 


-2, 


.2865587E+00 


6, 


.7419982E-02 


-1.3854520E-02 


1.6782558E-02 


137 


3, 


.4068802E+03 


2, 


.6582786E+01 


-2 


.2671584E+00 


9, 


.5942609E-03 


6.7850203E-03 


-9.2851399E-03 


136 


3, 


.4319311E+03 


2, 


.6778913E+01 


-2 


.3022740E+00 


2, 


.7985127E-02 


-8.5313181E-03 


-2.0850798E-03 


135 


3, 


,4573574E+03 


2, 


.6966587E+01 


-2 


.3222197E+00 


3, 


.5873646E-02 


-1.2983645E-02 


1.2479542E-03 


134 


3, 


.4831663E+03 


2, 


.7155877E+01 


-2, 


.3347319E+00 


3, 


.7873325E-02 


-6.2703785E-03 


1.1789081E-03 


133 


3, 


.5093650E+03 


2, 


.7334353E+01 


-2 


.3352488E+00 


1, 


.9924731E-02 


3.2293852E-03 


-7.4094784E-03 


132 


3, 


.5359529E+03 


2, 


.7538598E+01 


-2 


.3656308E+00 


3, 


.3747420E-02 


-9.2136909E-03 


-1.5741166E-03 


131 


3. 


.5629531E+03 


2, 


.7732403E+01 


-2 


.3769408E+00 


2, 


.9649622E-02 


-7.1513702E-03 


-1.9546620E-03 


130 


3, 


,5903666E+03 


2, 


.7933761E+01 


-2 


.3960270E+00 


3 


.4519896E-02 


-7.7112736E-03 


3.4520805E-04 


129 


3, 


.6182069E+03 


2, 


.8136593E+01 


-2 


.4112124E+00 


3, 


.4322472E-02 


-1.3197668E-02 


3.3426049E-04 


128 


3, 


.6464797E+03 


2, 


.8344177E+01 


-2 


.4316644E+00 


3, 


.6990870E-02 


-1.1417702E-02 


5.6447434E-04 


127 


3. 


.6751965E+03 


2, 


.8560555E+01 


-2 


.4575814E+00 


4, 


.3546365E-02 


-1.7266535E-02 


6.9978172E-03 


126 


3, 


,7043758E+03 


2, 


.8764023E+01 


-2 


.4596247E+00 


3, 


.2936489E-02 


-7.9666039E-03 


-7.1235356E-04 


125 


3, 


,7340186E+03 


2, 


.8980703E+01 


-2 


.4752705E+00 


3, 


.0423652E-02 


-8.3680162E-03 


-1.1760229E-03 


124 


3, 


.7641382E+03 


2, 


.9204603E+01 


-2 


.4968770E+00 


3, 


.4123134E-02 


-1.0928570E-02 


-7.7307545E-04 


123 


3, 


.7947486E+03 


2, 


.9427401E+01 


-2 


.5117336E+00 


3, 


.2840418E-02 


-7.2320179E-03 


-3.6479544E-03 


122 


3. 


.8258608E+03 


2, 


.9656143E+01 


-2 


.5307968E+00 


3, 


.1801252E-02 


-6.5872613E-03 


-2.3128226E-03 


121 


3, 


,8574873E+03 


2, 


.9886083E+01 


-2 


.5518506E+00 


2, 


.9777719E-02 


-8.6517922E-03 


-6.3210983E-03 


120 


3, 


.8896425E+03 


3, 


.0121318E+01 


-2 


.5652903E+00 


2, 


.9130129E-02 


-2.0933691E-03 


-6.0873269E-03 


119 


3, 


.9223364E+03 


3, 


.0368623E+01 


-2 


.5902555E+00 


3, 


.4140415E-02 


-1.3312468E-02 


-2.9620043E-03 


118 


3, 


.9555834E+03 


3, 


.0611776E+01 


-2 


.6114255E+00 


3, 


.6610362E-02 


-9.5989623E-03 


-1.7570264E-03 


117 


3. 


.9894010E+03 


3, 


.0858347E+01 


-2 


.6269035E+00 


3, 


.0456590E-02 


-3.8195133E-03 


-4.5113309E-03 


116 


4, 


.0237990E+03 


3, 


.1117735E+01 


-2 


.6542921E+00 


3, 


.6451959E-02 


-7.5083039E-03 


-2.7498142E-03 


Table 2 























The coefficients of the first 6 Legendre polynomials C0-C5 in the long-slit calibration 
lamp image for the 390 grating setting are listed in the table. 
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order 




Co 




Ci 




C2 




C3 


C4 


Cs 


139 


3, 


.7640555E+03 


2, 


.9247633E+01 


-2, 


.4983034E+00 


2, 


.1809631E-02 


-1.5302941E-03 


-7.3419534E-03 


138 


3, 


.7946624E+03 


2, 


.9476363E+01 


-2, 


.5209971E+00 


2, 


.7012675E-02 


-5.1708948E-03 


-3.3158614E-03 


137 


3, 


.8257730E+03 


2, 


.9708029E+01 


-2, 


.5427480E+00 


3, 


.1988803E-02 


-9.4145740E-03 


-1.3762068E-03 


136 


3, 


.8573989E+03 


2, 


.9940130E+01 


-2, 


.5606036E+00 


3, 


.3863295E-02 


-9.1725422E-03 


1.5817056E-03 


135 


3, 


,8895535E+03 


3, 


.0170043E+01 


-2, 


,5755437E+00 


2, 


.6638458E-02 


-7.7280717E-03 


-6.9431765E-03 


134 


3, 


.9222446E+03 


3, 


.0413128E+01 


-2, 


.5979035E+00 


3, 


.2158675E-02 


-9.7994044E-03 


-1.4903875E-03 


133 


3, 


.9554912E+03 


3, 


.0656432E+01 


-2, 


.6178296E+00 


3, 


.3148307E-02 


-6.8278925E-03 


-1.5417298E-03 


132 


3, 


.9893084E+03 


3 


.0903140E+01 


-2, 


.6327016E+00 


2, 


.7640170E-02 


-5.7173851E-03 


-4.5849886E-03 


131 


4. 


.0237052E+03 


3, 


.1159181E+01 


-2. 


.6591523E+00 


3, 


.4567971E-02 


-7.0897964E-03 


-2.5973397E-03 


130 


4, 


.0587017E+03 


3, 


.1416673E+01 


-2, 


,6778756E+00 


3, 


.5027143E-02 


-9.9222133E-03 


-2.2467367E-03 


129 


4, 


.0943163E+03 


3, 


.1668129E+01 


-2, 


.6874486E+00 


2, 


.3789638E-02 


9.4800526E-04 


-8.3427076E-03 


128 


4, 


.1305553E+03 


3, 


.1948742E+01 


-2, 


.7185913E+00 


4, 


.1422843E-02 


-8.0693347E-03 


2.1152305E-03 


127 


4. 


.1674438E+03 


3, 


.2213865E+01 


-2, 


.7383569E+00 


2, 


.9673746E-02 


-3.1165309E-03 


-7.7851206E-03 


126 


4, 


,2049961E+03 


3, 


.2492623E+01 


-2, 


,7655304E+00 


3, 


.4471338E-02 


-8.6897107E-03 


-5.5739705E-03 


125 


4, 


,2432354E+03 


3, 


.2768542E+01 


-2, 


,7788524E+00 


2, 


.7706485E-02 


-2.3526557E-03 


-5.2609012E-03 


124 


4, 


.2821701E+03 


3, 


.3066335E+01 


-2, 


.8149856E+00 


3, 


.9617907E-02 


-1.1644463E-02 


-6.7238977E-04 


123 


4, 


.3218291E+03 


3, 


.3357925E+01 


-2, 


.8360519E+00 


3, 


.8452345E-02 


-1.1566467E-02 


-1.5868857E-03 


122 


4, 


.3622303E+03 


3, 


.3658196E+01 


-2, 


.8606765E+00 


3, 


.7344307E-02 


-7.4652479E-03 


-5.7107517E-03 


121 


4, 


,4033926E+03 


3, 


.3960171E+01 


-2, 


,8889830E+00 


3, 


.1228066E-02 


-1.1341092E-02 


-7.0300975E-03 


120 


4, 


.4453403E+03 


3, 


.4275396E+01 


-2, 


.9108898E+00 


3, 


.7639151E-02 


-1.1469007E-02 


-1.2729088E-03 


119 


4, 


.4880943E+03 


3, 


.4594699E+01 


-2, 


.9380710E+00 


4, 


.0701831E-02 


-8.3199880E-03 


1.7192203E-03 


118 


4, 


.5316788E+03 


3 


.4918293E+01 


-2, 


.9638680E+00 


3, 


.8092495E-02 


-1.1269793E-02 


-3.9247857E-03 


117 


4. 


.5761167E+03 


3, 


.5253461E+01 


-2, 


.9941438E+00 


4, 


.5003528E-02 


-1.1020847E-02 


3.4696559E-03 


116 


4, 


.6214372E+03 


3, 


.5589023E+01 


-3, 


.0209908E+00 


3, 


.9108559E-02 


-6.1143803E-03 


-2.0235974E-03 


Tabic 3 























The coefficients of the first 6 Legendre polynomials C0-C5 in the long-slit calibration 
lamp image for the 430 grating setting are listed in the table. 
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Grating Setting 


Grating Setting 


390 


430 


Order Std. Dev. Largest Err. 


Num. Lines Order Std. Dev. Largest Error Num. Lines 



139 


5.40E-07 


9.54E-07 


13 


123 


1.16E-07 


2.57E-07 


16 


138 


3.72E-07 


-6.77E-07 


14 


122 


1.99E-07 


3.52E-07 


22 


137 


2.08E-07 


-3.20E-07 


14 


121 


2.61E-07 


4.94E-07 


19 


136 


3.62E-07 


7.22E-07 


17 


120 


2.35E-07 


-4.85E-07 


20 


135 


2.67E-07 


-4.98E-07 


21 


119 


2.35E-07 


-4.28E-07 


25 


134 


1.98E-07 


4.12E-07 


20 


118 


3.74E-07 


-5.75E-07 


18 


133 


1.65E-07 


3.83E-07 


15 


117 


1.65E-07 


3.49E-07 


21 


132 


1.06E-07 


-2.47E-07 


16 


116 


1.82E-07 


-4.30E-07 


22 


131 


1.83E-07 


-3.06E-07 


21 


115 


2.70E-07 


5.56E-07 


23 


130 


3.92E-07 


-7.24E-07 


24 


114 


1.58E-07 


-3.45E-07 


23 


129 


2.38E-07 


3.77E-07 


19 


113 


1.59E-07 


3.04E-07 


13 


128 


2.59E-07 


-4.93E-07 


15 


112 


3.29E-07 


-5.97E-07 


21 


127 


1.96E-07 


4.24E-07 


25 


111 


1.90E-07 


-3.97E-07 


23 


126 


2.78E-07 


4.72E-07 


24 


110 


2.44E-07 


-4.14E-07 


21 


125 


2.15E-07 


4.70E-07 


20 


109 


1.96E-07 


-3.13E-07 


20 


124 


2.64E-07 


4.44E-07 


25 


108 


1.99E-07 


-3.63E-07 


17 




Z.o 1 rLz-U 1 




zo 






-4i:.uyii/-u 1 




122 


2.81E-07 


-5.08E-07 


19 


106 


1.76E-07 


-2.86E-07 


16 


121 


2.42E-07 


4.41E-07 


19 










120 


2.56E-07 


-4.33E-07 


21 










119 


2.58E-07 


-4.51E-07 


19 










118 


2.26E-07 


3.26E-07 


18 










117 
Table 4 


1.72E-07 


-2.78E-07 


21 











Statistics for the fractional errors in wavelength for the 390 and 430 grating setting 
wavelength solutions. 
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UVES 



firstpassj 
First pass at order defin 


ongcal.pro 

tion, extractiDn and line 




waverefins 
Interactively sets good c 

Ma 


' uves.pro 

^1 lines and wavetength 
lion' 



-Seiasrid jlS^ on master long slit cali!bra8on:and 



,Ma£t£tr calibration line 
^rugture 



1 



Master data tiles for each grating angle 

Vi^vslength solution 
OalibrErtton line spectrum 
calibration line image 
line structure 
flat field correction image 
:gj*ter trace file 
.^«#umst)iflsfile 
slit yrfwid&i file 



Fig. 1. Overview flow chart for the wavelength calibration process. The names that 
end in .pro at the top of the boxes are the names of the IDL based procedures that 
were developed by the authors. 
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fi rstpa ss_la ng ca I - pro 



ASCII 1il8 of 
urns file 
names 



filesort.pro 

Softs file rajnes into otjservation categories 



mskemedian.pro 

Maite a mediar of ail long slit ra^bratioti lamp and OFder defirrtion images to nemoye cosmic lay hits 



(iefguess.pno 
Maizes a first pass at the order positions 



Output is 3 file of [orders .6|, [t>rd,0:ZJ are the 
i^Liadraiii; fit ccidclents, [citcl.j:^] sts the ^tan and 
and of ths ofder in pixals. 



ordtrace_uves.pro 
WskSs Ehesecofid and final determination of tiie enjSr positions 



Output is an array [nK.orrlers 2], 
[n!(.ori),0] = ypos. [n.i(,ortl,1] =fwhrTi 



makeflat_in'e3. pro 
1^/3 kes a flat itriEige frtim ibs medisn of the kxig slit flats 



straylgtit_Lves. pro 
Renrows strsy l^ht from ttte long slit cal lamp median 



spece)ctract_Lves. pro 
EKtracts tfie long slit cat iamp spectrum into indiuidual strips in each order 



Output it an array (if 
[rK,c*d,6.long)wiiere long is 
tt» Isngtti of ttie ejttf action . 



a£samtileltingcal_usrss.pro 
Connbines 'tis apectral strips into ons spectrum 



Wavelength solution 

derived from Die 
U\'£S pipeline data 
for fiach order 



Output is a spectrum of [nK,ofd,6j wtiere ttw 
firrst 3 elements are pijiel n urnber, flux , and 
noise, otlier 3 populated later, 





waveputlongcal.pro 
Put^ in Die initial guess of the wavelength Solution 


* 


Output Is a spatlrum of [n>L,o(i3,61 wiiere the 
fist 3 eiemnants are vacuunn wavelength, flux 
and noise, cltier 3 populated later 







flii£tlliis£tnjctijre_uvas ,ptio 
Producs a Ttcst guess line stajcture 




Output Is an IDL sauaut^ containing Una 
parameters and positions 






Fig. 2. Flow chart for the first pass analysis of the long-slit calibration lamp spectrum 
and wavelength calibration. The names that end in .pro at the top of the boxes are 
the names of the IDL based procedures that we developed. 
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Fig. 3. The order definition image for UVES at tlie 390 grating setting (centered at 
390 nm). Note the variation in intensity in both the horizontal and vertical directions 
and the slight upward tilt of the orders. 
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secon d pa s sj ong ca t . p ra 




I . 

Mal(e_lon3cal_nnask. pno 
Firtd the good cal lines 



Calcu^aie tie shifts fmm the cross conrslaticms 



Shift ttie images tif irtef polation 



Make the fmdian of the stuftad long cal images 



Output is an image of Sie 
iredian 



strayliglit._uir6s. pro 
Remme the sti^y light from the linage 



Output is ttw imstjer Idng silt calitiralion iamp 
image 



specs )iliact_uv6s.fm) 
Ertraois tire lorg slit cal ianp sjKctmm irto individuai slips ir each order 



Output is aiv airay of 
[rMijQiKj.S.loncf] whwSoii^ Is 
the lergth of ths extraction. 



assaml]jelongcai_uves. pro 
Combines the spectral strfps into one spectrum 



Output is a spectrum i3f[n».ord, 6] where the 
f rrst 3 elemer:ts at^ pixel number, II uy, and 
noise, other 3 populated later. 



waveputlorgoal. pn3 
ptils ir th& wavelerigth solirtion from (he diieetoiy 



Output ts a spactium of [ra.ond.S] whsrs tiiB 
ftrst 3 eiemerrts are vacuum wa^^la^gth, Ilux 
and noise, other 3 populated later 



fiiBtiineposition uvBS.pfo 




Oulpul is an IDL siriidture containing line 


Produce lha second line s.tructuw 




paran^ters ^nd positions 



lineJnfosKcahnga-pro 
Put In the line tise Inlarmalion found previously . 



Output Is an IDL structuns cxuntaining Una 
parameters, positions at>il use codes 



Fig. 4. Flow chart for the second pass analysis of the long-slit calibration lamp 
spectrum and wavelength calibration. The names that end in .pro at the top of the 
boxes are the names of the IDL based procedures that we developed. 
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1x10"® -5x10"^ 5x10"^ 1x10"® 

FRACTIONAL ERROR 



Fig. 5. The histogram of the hne wavelength fitting fractional residuals AA/A for 
the 390 grating setting 



Ld 

m 




5x10"^ 1x10"® 
FRACTIONAL ERROR 



Fig. 6. The histogram of the line wavelength fitting fractional residuals AA/A for 
the 430 grating setting 
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Fig. 7. The plot of all of the wavelength fitting fractional residuals A A/ A for the 
390 grating setting. The solid black line is the mean of the residuals in 100 A bins 
to check for systematic residuals. 
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Fig. 8. The plot of all of the wavelength fitting fractional residuals AA/A for the 
430 grating setting. The solid black line is the mean of the residuals in 100 A bins 

to check for systematic residuals. 

A A typical UVES pipeline reduction script used in this work 



The following is a set of typical commands is used with the MIDAS based 
UVES pipeline to produce the pipeline products used in this analysis. It is 
presented for completeness and as an example of the process. Individual batch 
reductions may vary slightly from the example. The commands are issued in 
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the directory that contains the data recovered from the archive. This script 
was for the Q0347-383 spectra at the 430 grating setting. 

fimidas 

confit/displ 

spht/uves refcalb.cat refcalb_spht.cat 
create/icat referb.cat null do -classification 
add/icat referb.cat thargood_2.tbl 
prepare/caldb refcalb_split.cat referb.cat 
split /uves object, cat object .split, cat 

reduce/uves object^plit.cat object j:edb. cat ref430.2x2.cat E optimal median 

After the script was run the appropriate bdf files were converted to FITS 
format using the outdisk/fits command. 
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